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SUMMARY 

A coordinated effort has been underway over the past four years to elevate unstructured -grid methodology 
to a mature level. The goal of this endeavor is to provide a validated capability to non-expert users for performing 
rapid aerodynamic analysis and design of complex configurations. The Euler component of the system is well 
developed, and is impacting a broad spectrum of engineering needs with capabilities such as rapid grid generation 
and inviscid flow analysis, inverse design, interactive boundary layers, and propulsion effects. Progress is also 
being made in the more tenuous Navier-Stokes component of the system. A robust grid generator is under develop- 
ment for constructing quality thin-layer tetrahedral grids, along with a companion Navier-Stokes flow solver. This 
paper presents an overview of this effort, along with a perspective on the present and future status of the methodol- 
ogy- 


INTRODUCTION 

The role of advanced Computational Fluid Dynamics (CFD) is becoming increasingly important in the 
aerodynamic design process. The transition of CFD technology from the research lab to the project environment 
has progressed rapidly in recent years with encouraging results. However, the intense competitive pressures on 
U.S. industries to build more efficient or specialized products, e.g. aircraft, automobiles, buildings, pumps, heat 
exchangers, etc., in less lime and at lower cost and risk continue to challenge the most ardent CFD method devel- 
oper. The challenge now is to complete the transition and realize the full potential of CFD technology in the engi- 
neering environment. The driving issues affecting CFD utilization in such an environment are time, confidence, 
cost, and risk [1]. Thus, codes to be applied in the project environment are those which are validated, robust, 
readily usable by non-CFD expert engineers, and provide reliable, timely information at a reasonable cost. 

The unstructured-grid (USG) approach has moved us much closer toward meeting these challenges. This 
methodology now enables modestly trained engineers to generate tetrahedral grids on complex configurations in a 
matter of days [2]. Inviscid flow solutions are readily obtainable from a number of good unstructured flow solvers 
[3]. Significant progress is also being made toward the grid generation and flow solution for viscous-flow problems 
[3], but additional work is needed before a reliable viscous USG capability emerges. This paper will include the 
status of our progress in this area. 

The status of USG Euler methodology is fairly mature and can play a significant role in the preliminary 
design phase of aerodynamic products. This inviscid technology has proven to be quite valuable as an engineering 
‘tool’ for providing critical insights which can lead to the timely solution of complex aerodynamic problems. For 
example, Ref. [4] describes a successful effort to redesign the wing/pylon intersection junction of a transport air- 
craft with under-wing nacelles by reducing the inviscid peak suction pressures in the region prone to viscous flow 
separation. The resulting drag reduction was verified by flight test. Work is also underway to assess the utility of 
this technology for modeling the underhood flow of automobiles. 
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The utility of USG Euler technology has been further enhanced by coupling it with other established engi- 
neering tools. As will be described in the text, an inexpensive 'viscous' capability is provided by coupling the flow 
solver with an interactive boundary-layer method. A design capability is available with the addition ot a powerful 
inverse design algorithm. Propulsion ctlecLs can be included through special boundary conditions. And while not 
discussed in this paper, work is underway to link a USG flow solver and grid generator with an acroelastic module. 
When complete, an engineering capability should be available for the rapid design of aerodynamic surtaccs for 
applicable complex configurations, which include the effects of propulsion flow, acroelastic deformation, and a 
boundary layer. 

There has been an ongoing coordinated effort at NASA LaRC over the past four years to mature the USG 
methodology. The goal is to provide a validated capability to non-expert users for performing rapid aerodynamic 
analysis and design ol complex configurations. This has been accomplished by advancing the methodology in a 
team environment using an application-oriented approach where the method development/calibration was guided 
by parallel focused application assessments from selected users within both NASA and industry. The remainder of 
this paper will present an overview ot our system, and a perspective on the present and luture status of the technol- 
ogy- 


OVERVIEW OF CURRENT CAPABILITY 

A diagram of the subject flow analysis system is presented in Figure 1. The primary components consist of 
a graphic user interface (GRIDTOOL), an advancing-front grid generator (VGRID), a cell-centered, upwind flow 
solver (USM3D), and a graphic post-processing code (VPLOT3D). 


Grid Generation 


GRIDTOOL. GR1DTOOL [5] is a graphic user interface (GUI) utility developed to serve as a “window” 

into the unstructured-grid generator, VGRID. The user must provide a surface definition in one of several formats, 
i.e. IGES, GRIDGEN, PLOT3D, LAWGS, or simple networks, which are converted within GRIDTOOL into a 
NURBS surface for manipulation. This GUI provides easy access to many functions, to be described below, which 
are useful in setting-up an input file for VGRID. 

VGRID . VGRID [2,6] is an interactive program for generating unstructured, tetrahedral/triangular grids 

by the advancing-front method [7]. A grid is generated by forming cells starting from the domain boundaries 
marching towards the interior ot the computational domain. The advancing-front process is illustrated in Figure 2 
for a simple 2D domain. The configuration ot interest is first defined in terms of a number of surface patches (line 
segments in 2D). A transparent ‘background grid' is used to control the distribution of points. The boundary seg- 
ments are then triangulated (divided) individually to form the surface mesh or initial front. Next, tetrahedral (trian- 
gular) cells are generated over faces on the front by introducing new points in the field. During this process, old 
faces are replaced by new ones (reshaping the front), and the front is advanced in the field until the whole region is 
filled with grid cells. 

For a 3D configuration, the geometry is divided into subdomains, using GRIDTOOL, on which 3- or 4- 
sided parametric surface patches are defined. The process of surface triangulation within VGRID involves trans- 
forming the 3D patch into a 2D parameter plane while preserving the overall projected shape of the 3D patch. After 
surface triangulation in 2D, the newly generated points arc then transformed back onto the 3D patch by interpolat- 
ing from its boundary curves. Depending on the patch topology, it is possible for points not to return exactly to the 
original surface after the transformation. As noted in Figure 1, the surface grid can optionally be passed through 
GRIDTOOL to project the points onto the surface. 

The distribution of grid points is controlled by a ‘structured' background grid [8], This transparent grid 
consists of Cartesian mesh overlaying the entire domain upon which the user prescribes point and line sources 
to impose the desired spacing distribution. Parameters are available to control cell size, and the direction and inten- 
sity of spatial variation. A smooth variation of spacing is achieved throughout the computational domain by solv- 
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ing an elliptic partial differential equation on the Cartesian mesh. The approach is analogous to modeling heat 
diffusion from discrete heat sources in a conducting medium. A smooth variation of cell spacing is critical to the 
robustness of the grid generator and the quality of the grid. 

The volume grid can be generated on a workstation through a series of restarts with VGRID. The code 
exploits a recurrent local/global renumbering system [2] which drastically reduces the computer memory require- 
ment. In each run, a prescribed number of grid points and cells are generated. With the renumbering system, the 
point and cell numbers, in each restart run, start from one rather than the accumulated global numbers. This allows 
specification of relatively small array dimensions regardless of the final grid size. 

After VGRID has completed the volume grid, there are usually some distorted tetrahedra, or some voids in 
the domain where the grid generator could not close the front. These problems occur when fronts approach from 
different directions in such a way that an ideal next point cannot be identified. As noted in Figure 1, the generated 
volume grid is postprocessed by the code POSTGRID [2] which removes distorted cells and/or a layer of cells 
around the problem areas, creates new fronts, and resumes the advancing-front process to close the grid. Since the 
modified front is now smoother after removing the anomalies, the problems are usually corrected, and a better 
quality tetrahedral mesh is produced. 

The new versions of VGRID and POSTGRID incorporate recent improvements for better grid quality, 
robustness, and easier grid generation. The enhancements include improved graphics for better visualization of 
interactive surface and volume grid generation, graphic version of POSTGRID, interactive check and improvement 
of the surface and volume grid quality, dynamic memory allocation, better CPU efficiency, and many other 
improvements. VGRID requires about 80 bytes of memory per tetrahedron, and generates about 220 tetrahedra 

(including I/O) per second using a Silicon Graphics Indigo 2 workstation with R4000 processor. 

Flow Solver 

USM3D . — USM3D [9-11] is a cell-centered, finite-volume upwind flow solver for solving the Euler equa- 
tions on tetrahedral grids. Inviscid flux quantities are computed across each cell face using the Roe [12] flux-differ- 
ence splitting approach (FDS), or the Van Leer [13] flux-vector splitting technique (FVS). Spatial discretization is 
accomplished by a novel cell reconstruction process, which is based on an analytical formulation for computing 
solution gradients within tetrahedral cells. Solutions are advanced in time by either a 3-stage Jameson-style Runge- 
Kutta explicit time-stepping scheme with convergence accelerated to steady state by local time stepping and 
implicit residual smoothing, or the linearized backward Euler implicit scheme of Ref. [14]. 

Cell reconstruction scheme : The higher-order reconstruction scheme, derived in Ref. [11,15] and illus- 
trated in Figure 3, is based on a Taylor series expansion of the cell-averaged solution to the cell face. A key compo- 
nent of the scheme is the reconstruction of surrounding cell-averaged data to a common vertex or node by a 
weighted averaging procedure. Reference [9] proposed a scheme based on an inverse-distance weighted averaging 
of the primitive variables from the cell centroid to the cell vertices. While this approach has proven to be both 
accurate and robust through wide application, it is not fully second-order accurate in space. It has been shown in 
Ref. [16] to be approximately 1.85-order accurate. 

As development efforts progressed toward solving the Navier-Stokes equations on highly stretched tetra- 
hedral grids, it became evident that the accuracy of the inverse-distance averaging scheme was not adequate. Thus, 
a fully second-order accurate averaging procedure was implemented which is based on work by Holmes and Con- 
nell [17] and Rausch, et. al. [18]. The procedure is derived by solving a constrained minimization problem to deter- 
mine weight factors which satisfy Laplacian relationships. The algorithm reconstructs to machine accuracy the 
exact values of a linear function at a node from surrounding cell-centered function values on an arbitrary tetrahe- 
dral grid. Furthermore, the simple universal formula shown in Figure 3 for expanding the cell-centered data to the 
cell faces also reconstructs the exact value of a linear function to the cell face. Thus, the entire spatial reconstruc- 
tion scheme is termed second-order accurate, which has been verified by Mitchell [16]. 
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Limiter device : Euler, finite-volume, upwind methods occasionally encounter circumstances where some 
local cell pressures, p ce u y become negative and preclude convergence of the global solution. This often occurs in 
regions where the inviscid assumptions are violated and viscosity would normally introduce enough diffusion to 
relieve the problem in real flow. Since the Euler equations are inviscid, the only source of diffusion in the solution 
comes from the discretization error inherent in the numerical approximation. A special empirical technique has 
been developed to resolve these local anomalies as they occur and permit the global solution to advance toward 
convergence. During each time iteration cycle, p cen is not permitted to become less than a user prescribed mini- 
mum value, p toh i*. p cell =max(p celh p tol ). Furthermore, for higher-order differencing, cells which satisfy p ce u<p to i 
are set to and maintained with a more diffusive lst-order stencil by imposing <J> ce //= 0, as shown in Figure 3. In 
many cases, this technique introduces enough local numerical diffusion to stabilize the global solution and permit 
convergence. It has greatly increased the robustness of USM3D for complex geometries. 

Data structure : An efficient face -based data structure [19] is exploited in the flow code, which does not 
require interface fluxes to be stored during the solution cycle. The flux algorithm cycles through the triangular 
faces and accumulates the residuals by scattering each flux to the two adjacent tetrahedral cells. For proper vector- 
ization, face-based data structure requires that the faces be ordered into groups or “colors” in which no two faces 
share the same cell or opposing node. The explicit time integration scheme requires 45 words/cell of memory and 
18 ps/cell/cycle (higher order) on a Cray Y-MP supercomputer. Similarly, the implicit scheme requires 180 words/ 
cell and 64 ps/cell/cycle (higher order) on the Cray Y-MP. While the implicit scheme requires more CPU time per 
iteration than the explicit scheme, it is demonstrated in Ref. [11] that the implict algorithm converges in 6-times 
less computer time than the explict one since higher Courant-Friedrichs-Lewy (CFL) numbers can be utilized. The 
following table shows a machine-to-machine comparison of Central Processing Units (CPU) requirements for 
USM3D on a test grid for the ONERA M6 wing. The computations involved a combination of 1st- and higher- 
order differencing. No special effort was made to enhance vectorization or performance on the non-Cray machines. 

The larger relative CPU time for the Iris Indigo 2 between implicit (98.9) and explicit (27.6) time integration is 
related to excessive swapping of the buffer due to insufficient core memory. Implicit calculations were not made on 
the Sun 4 because of insufficient memory. 
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Boundary Conditions . — The USM3D code is designed for the easy addition/modification of boundary con- 
ditions (B.C.). It contains standard B.C.’s for flow tangency or no-slip on solid surfaces, characteristic inflow/out- 
flow for subsonic outer boundaries, and freestream inflow and extrapolation outflow for supersonic flow. Some 
additional special boundary conditions, which have been added to improve the robustness of the code or to expand 
its capability, are described in the following. 

Blunt base : With the application of USM3D to a wide range of complex 3D configurations, numerical dif- 
ficulties routinely occur in blunt-based regions where a wake-type flow would develop in real flow, e.g. a blunt 
wing trailing edge or a blunt termination of a fuselage. For such cases, the inviscid flow assumption yields a singu- 
larity where the flow attempts to negotiate the corner. 
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A special boundary condition has been developed to mimic the relieving effect of a blunt-base wake on the 
surrounding inviscid flow field, which is to provide a smooth departure of the flow past the corner. As illustrated in 
Figure 4, the reconstructed primitive variables, (pH,v,w,/>), from the nodes of blunt-base boundary cells are aver- 
aged and assigned directly to the boundary face and ghost cell. The effect is to introduce a solution defined transpi- 
ration boundary condition to the subject faces which simulates the presence of a wake. 

Engine model : USM3D contains a feature for modeling propulsion-induced effects of up to four jet-engine 
or jet-nozzle components. The model [20] computes inflow and/or outflow boundary conditions from user pre- 
scribed values for the total pressure and temperature change across or within the combustor, the exit pressure, the 
fuel fraction, and the direction cosines for the jet. When applicable, the inlet mass flow is automatically balanced 
with that of the nozzle. 

Preprocessing. — A preprocessing step is involved before running USM3D on a given grid generated by 
VGRID. The VGRID code outputs a coordinate file of the vertex locations, a cell-to-node connectivity file, and a 
triangle-to-patch correlation file. The PREFLO code, depicted in Figure 1 generates an additional face-to-cell con- 
nectivity file which is required for the face-based data structure in USM3D. (There are plans to modify the VGRID 
output to generate the appropriate files and eliminate this preprocessing step.) Boundary conditions are applied at 
the time of USM3D execution through a B.C. map file output by GRIDTOOL. 

Graphic Postprocessing, VPLOT3D 

VPLOT3D is an interactive, menu-driven post processing program for extraction and display of data on 
unstructured (tetrahedral) grids. The philosophy behind the development of VPLOT3D has been to keep the pro- 
gram as simple and straight-forward as possible. There is only one ‘text’ window and only one main ‘graphics dis- 
play’ window. Most of the interaction and steering takes place through pop-up menus. 

VPLOT3D has the ability to display grids and flow quantities on either boundary surfaces or user defined 
arbitrary planes in the field. Flow quantities can be displayed using line or filled contours, velocity vectors and par- 
ticle traces. Some unique features of the program include dynamic memory allocation, a ‘probe’ option that lets a 
user query the local value of a certain displayed fluid dynamic quantity, ability to plot two different data sets side- 
by-side and a journaling (log file) capability. Grid post-processing options allow a user to isolate and display ‘bad’ 
grid cells and to visually confirm boundary conditions before proceeding for flow solution. Several hardcopy out- 
put options are available including Black&White and Grayscale Postscript and RGB formats. 

TECHNOLOGY ASSESSMENT 

The present USG flow analysis system is being widely used throughout industry and government laborato- 
ries. Studies have been published related to subsonic transports [4,21], high-speed research [22,23], military air- 
craft technology [24], turbomachinery applications [25,26], and code calibration/assessments [27-29]. Results 
from a published code calibration study are included herein to illustrate the capability. 

A study was performed to assess the accuracy and utility of VGRID/USM3D on a wing/pylon/finned store 
configuration [27] shown in Figure 5. Tetrahedral grids were generated for the configuration with baseline (or 
store-in-carriage) store position, two dropped positions for M^-0.95, and two positions for M w =1.2. The initial 
grid was generated prior to the development ot GRIDTOOL and took 4 days to complete. (The same task would 
require 1 or 2 days with present codes.) The remaining four grids were generated by ‘freezing’ the baseline surface 
triangulation while translating/rotating the surface patches, background-grid sources, and surface grid for the store 
to prescribed locations, then generating a new volume grid. The latter four grids required approximately one-hour 
each to generate. 

Solutions were computed with Runge-Kulta time stepping and convergence acceleration using a CFL 
number of 4 in about 2000 iterations, while the L2 norm (rms average of all residuals) decreased by about 4 orders 
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of magnitude. For the baseline position, the computation required about 5.7h of Cray Y-MP time, whereas for the 
other two positions, the required computational time was of the order of 4 CPU h. 

Figure 5 compares computed and experimental Cp data [30] on the store along two streamwise rays (along 

top, <t>=5°, and bottom, $=185°) from the nose to the tail. For each case, the computed pressure coefficient com- 
pares very well with experimental data. In particular, note the comparison along the 5-deg ray for the baseline con- 
figuration that passes through the gap between the store and the pylon and where a large interference is expected. 
While not shown, very good agreement was also obtained on the fins of the store. 


STATUS OF EMERGING CAPABILITY 


Thin-Layer ‘Viscous’ Grid Generator 

Advancinp-Lavers Method . — A new method for generation of thin-layer ‘viscous’ grids has been recently 
introduced and successfully applied to 2D [31] and 3D [32] problems. The method, referred to as the advancing- 
layers method (ALM), is based entirely on a modified advancing-front method (AFM) and benefits from the gener- 
ality and flexibility of the conventional advancing-front-based Euler grid generators. Being based on an unstruc- 
tured grid strategy, the method alleviates the structural limitations of the structured and semi-unstructured 
techniques. 

The generation of a grid with the present approach is divided into three separate stages: 1) surface grid 
generation, 2) construction of high-aspect-ratio cells in the viscous region, and 3) generation of regular (isotropic) 
cells in the inviscid-flow region. Steps 1 and 3 utilize established methodology encompassed in the advancing- 
front inviscid grid generation code VGRID. The second step is performed independently by the ALM to reduce the 
overall complexity of the process. Although the generation of grid is divided into separate stages with the present 
approach, the entire process is performed in a single run with automatic transition from one stage to another. 

The main features of the advancing-layers technique are similar to those of the conventional AFM. A vol- 
ume grid is generated through a marching process in which tetrahedral cells originate from a triangular surface 
mesh and proceed into the computational domain. In contrast to the conventional method which adds cells in the 
field in no systematic sequence, the ALM advances one layer of cells at a time to reduce the complexity of produc- 
ing high aspect ratio cells. Grid cells are formed by connecting new points, inserted along predetermined surface 
vectors using a prescribed stretching function, to the corresponding faces on the front. During the process, the 
integrity of the generated grid is monitored at each step. The layers continue to advance in the field, while growing 
in thickness, until either 1) opposite fronts approach to within a local cell size or 2) certain grid quality criteria, dic- 
tated by a global background grid, are locally satisfied. When the conditions are met on all faces on the front, the 
process automatically switches from the advancing-layers to the advancing-front mode to generate equilateral cells 
in the remaining inviscid-flow region. 

To demonstrate the capability of the method, sample 2D and 3D grids around complex multi-component 
configurations are presented. The 2D geometry consists of three airfoil sections in a high-lift arrangement. The 
generated grid, shown in Figure 6, contains 34,987 cells and 17,705 points with a first layer spacing of about 

7X10* 6 main airfoil chord length. The flow solution, shown in Figure 6 was obtained on the grid at M^-0.2, 
a=16.2°, and Re^= 9X10 6 using an available node-based, upwind flow solver [33] with the Baldwin-Barth turbu- 
lence model. 

A representative complex 3D configuration was constructed to show the ability of the method to negotiate 
difficult aspects of 3D viscous grid generation. Although the geometry resembles a quasi-two-dimensional shape, it 
contains 3D complexities such as multiple components, sharp convex and concave edges, small gaps, and thin sur- 
faces at the wing tip. The generated grid (surface and volume) is shown in Figure 7. The ‘viscous’ portion of the 
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grid contains 490,266 tetrahedral cells and has been generated using an SGI Indigo 2 workstation (R4000 proces- 
sor) in only 4.3 minutes. The final grid contains a total of 168,734 points and 958,716 cells. 

Anisotropic stretching: — In order to generate viscous grids of reasonable size on complex configurations, 
it will be necessary to have anisotropic (stretched) surface grids so that fewer points can be distributed in directions 
of reduced flow field gradients. The development of such a capability is currently under way which has produced 
satisfactory preliminary results. Figure 8 compares two surface grids for a wing configuration: one isotropic and 
another stretched along the leading and trailing edges of the wing. Also the field grid distributions are represented 
by triangulating a plane passing through the middle of the wing. The isotropic surface mesh (wing only) contains 
40,778 equilateral triangles and is shown in Figure 8(a). A corresponding stretched grid with a similar chordwise 
grid distribution is shown in Figure 8(b). As indicated, the grid is stretched along the edges of the wing in different 
directions, with the stretching propagated in the field close to the body and smoothly dissipated further away from 
the geometry. The maximum cell aspect ratio in this grid is about 20:1 at the leading edge of the wing. Although 
the stretched grid has a similar distribution of points in the chordwise direction, it contains only 6,841 triangles on 
the wing. 


Navier-Stokes Flow Solver 

Considerable progress has been made toward solving the Navier-Stokes equations on thin-layer tetrahedral 
grids with USM3D. Ref. [11] provided an initial validation of the solution algorithm on a tetrehedral grid con- 
structed by subdividing the hexahedral cells of a structured grid. The following presents a full demonstration of 
computing a viscous-flow solution on a fully unstructured thin-layer tetrahedral grid produced from VGRID. 

Hummel delta wing . — An unstructured viscous grid of 730,454 tetrahedral cells was generated for a delta 
wing of aspect-ratio 1 by an ALM version of the VGRID code. Normal spacing was distributed with geometric 
stretching near the surface to yield approximately 10 nodes (30 tetrahedra) in the mid-chord boundary layer for 
laminar flow. Anisotropic stretching, evident in the surface grid shown in the Figure 9, was utilized to reduce the 
chordwise density of cells in the grid — this feature resulted in a 7-fold reduction in grid size compared to standard 
isotropic triangles. 

A viscous flow solution was computed with USM3D. A companion structured-grid computation was also 
performed with the CFL3D code [34] on a 65X65X33 (span-radial-chordwise) H-O grid for comparison. Both the 
unstructured and structured computations were performed on the delta wing at M 00 =^.3, a-20.5 ° , and 
Re L =0. 9X1 0 6 with laminar flow. 

The figure portrays the surface “oil-flow” pattern and a comparison of the spanwise distribution of pressure 
coefficient, C p , at four chord stations. The coalescing and diverging streamlines in the computed flow patterns 
show evidence of the primary, secondary, and tertiary vortices. Good agreement is noted between the unstructured 
and structured solutions, and with the experimental data of Hummel [35]. 

Turbulence modeling . — The one equation turbulence model of Spalart-Allmaras [36] has been installed in 
USM3D and is presently undergoing testing. The derivation of this model relies heavily on empirical reasoning 
and dimensional analysis, and is receiving a growing acceptance within the CFD community. The turbulence equa- 
tion solves for turbulent viscosity separately from the flow equations during each time step, resulting in a loosely 
coupled solution process. The governing equation is solved using the same backward-Euler time-stepping scheme 
[14] as used for the flow equations. 

A preliminary assessment of the turbulence model has been made for flat-plate boundary layer flow. A 
two-dimensional flow was simulated by computing on a three-dimensional ‘channel’ grid which was constructed 

by subdividing the hexahedral cells of a structured grid into tetrahedra. The grid was sized for a y + « 1 and 20-25 
nodes across the boundary layer with M^~0.5 and Re L =2A70 6 . Figure 10 depicts a very good agreement between 
the computed behavior of the law-of-the-wall and the model of Spalding [37]. Additional assessments are under- 
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way for a NACA 0012 airfoil and the ONERA M6 wing. The use of wall functions to eliminate the need to resolve 
the laminar sublayer is also under investigation. 

Interactive Boundary Layer 

An interactive boundary layer (IBL) capability has been installed into USM3D [38] to provide for an inex- 
pensive accounting of viscous effects in an unstructured-grid flow solver for predominantly attached flow cases, 
such as the transonic cruise condition of a subsonic transport. The two-dimensional strip boundary-layer method of 
Cebeci et al. [39] for turbulent, compressible flow was integrated with USM3D in a coupled fashion where the dis- 
placement thickness was simulated by transpiration velocities. The automated process begins with a partially con- 
verged inviscid solution and continues with the boundary layer method being invoked at specified intervals 
throughout the remainder of the solution iterations. The inviscid pressures along chordwise strips on the wing sur- 
face are passed to the boundary layer routine which calculates the estimated boundary layer displacement thick- 
ness. The transpiration velocities are computed along each strip from the boundary layer displacement thickness 
using Lighthill’s equation [40], A linear interpolation is performed between the strips to establish transpiration 
velocities at the remaining surface nodes. The face-based displacement boundary conditions are determined by 
averaging the node-based transpiration velocities to the centroid of the boundary triangle, then modifying the orig- 
inal tangent-flow condition to include the non-zero normal velocity. Face-center density and energy are set equal to 
their corresponding cell-center value. A similar capability tor a node-centered unstructured flow solver is described 
in Ref. [41]. 

An application of the unstructured IBL capability, taken from Ref. [38], is shown in Figure 11 for a low- 
wing transport configuration with a supercritical airfoil at M^O.77, a=1.12°, and Re~=2.6Xl(ft. The computa- 
tional grid contains 286,000 tetrahedral cells and has ten chordwise strips along the wing span for applying the 
boundary layer method. A comparison of pressure coefficient at the 34-percent semispan station is shown in Figure 
11. The inclusion of the boundary layer model results in substantial improvement in predicting both shock location 
and strength, as well as in the suction plateau. The poor agreement between the Euler result and the experimental 
data is typical of inviscid solutions on supercritical airfoils, where the shock location can be overpredicted by as 
much as 10- to 20-percent aft-chord. 


Inverse Iterative Design 

The Constrained Direct Iterative Surface Curvature (CDISC) design method [42] has been coupled with 
USM3D to provide a capability for aerodynamic design of complex configurations. The CDISC method automati- 
cally generates and modifies the target pressures that are used by the basic DISC design module so that flow and 
geometry constraints are met. The CDISC module is coupled with the flow solver in an iterative manner and is 
called after a user-specified number of flow iterations (typically 25). After the target pressures are altered to meet 
the design constraints, the basic DISC module is used to modify the affected surface points. Geometry constraints 
are directly enforced on this new surface, then the volume grid is perturbed, using the spring analogy of Ref. [43], 
to reflect the geometry changes. Since the target pressure modification, surface design, and flow solution are con- 
verged in parallel, and the CDISC design system components themselves are very efficient, a converged design is 
usually obtained in about 25-50% more time beyond the initial baseline analysis. 

To illustrate the USM3D/CDISC design capability, the wing of an executive transport configuration (Fig- 
ure 12) has been redesigned for reduced drag at a primary design point without increasing the drag at a secondary 
design point. The baseline configuration was first analyzed using the USM3D code to establish reference values for 
performance comparisons at the primary (Mach = 0.8) and secondary (Mach = 0.9) design points. These conditions 
correspond to the long-range and high-speed cruise points, respectively. Two approaches are available in CDISC 
for reducing the wave drag of a configuration at transonic speeds: 1) a simple incremental method that limits the 
Mach number ahead of a shock to a specified value while matching lift and pitching-moment constraints, and 2) 
the Modified Uniform Distribution (MUD) method, which redistributes the chordwise loading over an airfoil sec- 
tion. The MUD option was selected for this case and was used to redesign the wing to reduce the drag at the pri- 
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mary design point. The wing planform, nacelle location, root airfoil section, airfoil maximum thickness, and 
spanwise load distribution were all held fixed. Initial results showed a 15 count drag improvement at Mach = 0.8, 
with only a 3 count penalty at Mach = 0.9. This design however had an unacceptable twist distribution just out- 
board of the nacelle. An input parameter to the MUD method was adjusted to shift the chordwise loading forward 
at several sections, which brought the twist distribution much closer to that of the original wing. The baseline and 
modified pressure distributions at Mach = 0.8, along with the corresponding airfoil sections, are shown in Figure 
12 for two wing stations. At the inboard station, the strong shock near the quarter-chord of the baseline wing has 
been significantly weakened, resulting in an overall drag decrease of 6 counts. Analysis at the high-speed cruise 
point indicated that the new design had a 1 count drag improvement over the baseline at that condition. The design 
took an additional 4 hours of Cray C-90 time beyond the 16 hours required to obtain the initial baseline analysis. 
The grid had 700,000 tetrahedral cells, producing a central memory requirement of 148 megawords. 

Solution Adaptive Gridding 

One of the strong advantages of unstructured-grid technology is the potential to exploit its flexibility in 
adapting or clustering the grid in the vicinity of dominant flow features, such as shocks and vortices, to enhance 
solution accuracy with fewer cells. This capability should facilitate a more rapid transfer of the technology to the 
workstation environment. Attention has been focused on providing an unstructured “remeshing” capability which 
distributes grid points efficiently throughout the domain and adequately resolves important flow features. 

Remeshing procedures have been developed and applied for shock [44] and vortex [45] dominated flows. 
The procedure integrates the grid generation and flow solver functions in such a way that the information produced 
by one is successively used by the other until a desired accuracy is obtained. Accordingly, the grid generator 
VGRID, is combined with the flow solver USM3D, to produce a unified procedure that produces input parameters 
for the mesh generator based on the solution on a previous grid. The concept of a background grid, which is used in 
the advancing-front technique to specify and control grid size distribution, is exploited in the remeshing procedure. 

Shocked flows : — The procedure of Ref. [44], which addressed shock dominated flows, was based on the 

“tetrahedral” background grid 1 concept [7] where the computational grid serves as the background grid for gener- 
ating a subsequent grid. Grid spacing parameters were calculated at the nodes of this mesh based on the previous 
flow solution, and a new mesh was generated. Figure 13 shows the result of applying this remeshing procedure on 
the ONERA M6 wing. The upper portion of the figure compares the upper surface triangulation between the 
adapted grid (at the end of three adaptation cycles) and an unadapted grid. On the lower portion, C p distribution at 

two spanwise stations are compared (for M oo =0.84 and a=3.06°) between the adapted and unadapted grids, and the 
experimental data. While both sets of computed data compare well with the experimental values, the efficiency of 
the adaptive procedure is established by the fact that the adapted grid had about 2.5 times fewer cells than the 
unadapted case. 

Vortex flows : — A solution adaptive remeshing procedure, based on the “structured” Cartesian background 
grid concept, has been developed for low-speed vortex dominated flows [45], such as those produced by sharp 
swept-back leading edges or shed from the free ends of lifting surfaces. As discussed earlier with the structured 
background grid, grid spacing is controlled by point and line ‘sources’ strategically placed throughout the domain 
by the user. The present approach provides a capability for automatically prescribing line sources for refining the 
grid in regions of vortex flow. 

In the current procedure, the vortices are located by an entropy parameter. In addition, a vortex in close 
proximity to a surface induces a static pressure peak on the surface. To accurately resolve both the vortex core and 
the accompanying static pressure peak, line sources are placed in the field and on the surface of the configuration. 
The procedure is automated to the point that the user need only specify the approximate starting locations for the 


1. The “tetrahedral” background grid is no longer used in VGRID, but has been replaced by the “structured” background grid 
approach. 
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vortices and threshold cut-off values for both the entropy and the static pressure. All other functions, such as iden- 
tifying the location and strength of the background grid ‘sources’ and creation of a new input for the grid generator, 
are automated. 

The adaptation procedure has been applied to the Modular Transonic Vortex Interaction (MTVI) model, 
tested at NASA Langley, which is characterized by the presence of multiple vortices (see Figure 14). Computations 
were initiated on a relatively coarse grid with 35,388 points and 188,304 tetrahedral cells (GRID 1). Two succes- 
sively finer adapted grids were generated with the adaptive procedure. The final adapted grid (GRID 3) had 68,158 
points and 371,360 cells. 

The left side of Figure 14 shows a comparison of surface triangulation for the unadaptcd (GRID 1) and the 
adapted (GRID 3) cases. For the adapted case, surface grid clustering near leading edges as well as near the wing 
tip is evident representing the static pressure peak due to proximity of vortices to the configuration surface. The 
bottom right portion of the figure shows a comparison of the field grid projected on to a cross sectional plane at 93- 
percent of root chord. Grid clustering around the vortex cores between the fuselage and the vertical tail and near 
the wing tip can be seen clearly. On the top right of the figure, flow field results are compared between the 
unadapted (GRID 1), the adapted (GRID 3) and experimental data. The C p comparison is shown at a streamwise 

station located at 93-percent of the root chord for M^-0.4 and a=10.54°. Results are also compared from a ‘fine’ 
grid with 148,285 points and 825,469 tetrahedra . During generation of the ‘fine’ grid, no special effort was made 
to cluster grid points around vortex core locations. An expected higher suction peak resulting from grid adaptation 
is evident. The fine grid did not capture the suction peak at the station shown due to lack of grid resolution there, 
while the adapted grid automatically provided the grid resolution needed. This case clearly establishes the ability 
of the adaptation procedure to automatically provide adequate grid resolution where needed. The whole adaptation 
process required about 2.5 hours of CPU time on a Cray-YMP computer, compared to about 4.0 hours for the ‘fine 
grid. The adapted results show a better accuracy with about 2.2 times less number of cells and about 38-percent 
CPU saving compared to the ‘fine’ grid case. 

CONCLUDING REMARKS 

An overview has been presented of a maturing capability for performing rapid aerodynamic analysis and 
design of complex configurations using unstructured-grid methodology. Development priorities have been dictated 
by program needs from within NASA as well as those of industry partners. This technology also offers strong 
potential for impacting other non-aerospace disciplines, such as automobiles, boat sail and hull design, wind loads 
on buildings, design of ventilation systems, analysis of pumps and pipe systems, design of heat exchangers, and 
electromagnetics. The present unstructured Euler technology can be readily applied by modestly trained engineers 
to solve complex aerodynamic problems. Significant progress has also been made towards developing an unstruc- 
tured Navier-Stokes capability for complex geometries. 

Considerable utility can be derived from the base Euler technology. Often, the inviscid flow equations can 
provide the critical aerodynamic insights needed to resolve more complex problems. The Euler methodology has 
also been enhanced with additional engineering features such as interacting boundary layers, constrained inverse 
design, propulsion flow boundary conditions, and solution-adaptive gridding, thereby providing a powerful capa- 
bility for aerodynamic design and analysis to the project-level engineer. 

The capability for solving the Navier-Stokes equations on complex configurations using unstructured-grid 
methodology is still evolving. Such a flow solver has been developed and is presently being tested. The primary 
obstacle to a practical unstructured Navier-Stokes capability has been the lack of a robust thin-layer ‘viscous’ tetra- 
hedral grid generator. As demonstrated in this paper, the Advancing-Layers method shows great promise toward 
providing that capability. 
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Figure 1. - Flowchart of unstructured-grid flow analysis system. 
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Figure 2. - Advancing-front process for a simple 2D domain 
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Figure 3. - Higher-order stencil for tetrahedral cell-centered scheme. 
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Figure 4. - Schematic of blunt-base boundary condition 
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Figure 5. - - Surface triangulation on a wing/pylon/store configuration, and surface pressure comparisons on 
store along streamwise rays, A/ x =(9. 95. 
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Figure 6. - Unstructured triangular viscous grid and viscous-flow solution for multi-element airfoil atA/ oo =0.2, 
a=16.2°, and Re c =9Xl(ft: a) surface pressure distributions and b) velocity profdes. 



Figure 7. - Complete tetrahedral viscous grid on a multi-element wing shown on a plane. 
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(b) Multidirectional anisotropically stretched surface and field grid (6,841 triangles on wing). 


Figure 8. - Example of multidirectional anisotropic grid stretching for cranked-wing configuration. 
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’igure 9. - Laminar flow over aspect-ratio 1 delta wing for M ao =0.3, a=20.5°, Rc\ =0 ,9X 1 (/’ . Unstructured 
surface grid and “oil-flow” pattern (left), pressure profiles at four chord stations (right). 
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Figure 10. - Comparison of computed law-of-the-wall with analytical model for flat plate boundary layer flow 
M^O.5, Re L =2X10 6 . 






Figure 11.- Effect of interacting boundary layer on pressure predictions for a low-wing transport. A/^^O.77, 
a=1.12 0 , and Re c =2.6X10 6 
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Figure 12. - Application of USM3D/CD1SC design method to wing redesign of an executive transport 
configuration, M^O.80. 
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